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We study nonlinear interactions in a strongly driven optomechanical cavity, in regimes where 
the interactions give rise to resonant scattering between optomechanical polaritons and are thus 
strongly enhanced. We use a Keldysh formulation and self-consistent perturbation theory, allowing 
us to include self energy diagrams at all orders in the interaction. Our main focus is understanding 
how non-equilibrium effects are modified by the polariton interactions, in particular the generation 
of non-zero effective polariton temperatures from vacuum fluctuations (both in the incident cavity 
drive and in the mechanical dissipation). We discuss how these effects could be observed in the 
output spectrum of the cavity. Our work also provides a technical toolkit that will be useful for 
studies of more complex optomechanical systems. 
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I. INTRODUCTION 

The rapidly growing field of cavity optomechanics 
seeks to understand the interaction between photons 
and mechanical motion in driven electromagnetic cav¬ 
ities, hopefully in truly quantum regime [lj [2]. The 
past few years have seen many breakthroughs, including 
laser cooling of mechanical motion to the ground state 
0 0: and the optomechanical generation of squeezed 
light leaving the cavity m- Although the basic radia¬ 
tion pressure interaction between a mechanical resonator 
and cavity photons is intrinsically nonlinear, almost all 
the remarkable achievements in the field to date rely on 
working in strongly driven regimes where the dynamics 
is essentially linear. To see truly nonlinear effects in the 
simplest setting, one needs to achieve a single photon, 
single phonon optomechanical coupling g which exceeds 
both the mechanical frequency u>m and the cavity damp¬ 
ing rate n [SHU. With the exception of experiments 
using cold atoms mug this parameter regime remains 
challenging for experiments. 

Nonlinear effects which only require g ~ k can be 
achieved in slightly more complex settings where the non¬ 
linear interaction becomes resonant. This can occur in 
an optomechanical setup with two optical modes mug. 
Alternatively, this can occur in a standard single mode 
optomechanical cavity which is driven, such that the in¬ 
teraction becomes resonant in a basis of dressed states 
(so-called optomechanical polaritons). This occurs both 
in regimes of weak driving mm and strong driving m 
similar physics can also occur in membrane-in-the-middle 
style optomechanical systems [H]- Note that such strong 
driving regimes, where the drive-enhanced optomechani¬ 
cal coupling exceeds dissipative rates, has been achieved 
in several experiments mm- 

In this paper, we both expand upon and extend the re¬ 
sults previously reported in Ref. \T7, We consider a stan¬ 
dard single-cavity optomechanical system that is driven 
(possibly strongly), in a regime where the nonlinear in¬ 
teraction is resonantly enhanced. We again address this 



FIG. 1: (Color online) (a) A generic optomechanical cavity: 
a cavity mode of frequency u)c and damping rate k is driven 
by a laser at frequency <j.l, and is coupled via radiation pres¬ 
sure to a damped mechanical mode (resonance frequency ujm, 
damping rate 7 ). The temperature of the mechanical bath is 
T m . (b) The single-photon optomechanical interaction g can 
be made resonant in the basis of optomechanical polaritons; 
the resulting enhanced interaction shown schematically and 
described in Eq. | 8 |. 


system using a Keldysh formulation [221 [23: and per¬ 
turbation theory; we begin this paper by providing a 
more complete discussion of this approach to optome¬ 
chanics, emphasizing the subtleties involved in treating 
dissipation. Our work thus contributes to a growing body 
of work using the Keldysh technique to address various 
quantum optics contexts [MH25] . Unlike our previous 
study, we go beyond a simple lowest-order treatment and 
introduce a self-consistent perturbation theory. This cor¬ 
responds to an infinite partial resummation of self-energy 
diagrams, and extends the range of couplings and tem¬ 
peratures that we can address. 

Our work also address a new set of physical phe¬ 
nomena. Ref. H7| focused on understanding the cav¬ 
ity density of states (DOS), the quantity probed in 
optomechanically-induced transparency (OMIT) exper¬ 
iments m l29H3Tj . Here, we instead focus on the non¬ 
equilibrium state of this system: what are the effective 
temperatures of the optomechanical polaritons? Of par- 
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ticular interest is how vacuum fluctuations can result in 
non-zero effective temperatures. Such effects are often 
termed “quantum heating” [52H5P] ; the simplest exam¬ 
ple is the amplification of zero-point fluctuations by a 
parametric amplifier. Though often described using dis¬ 
parate terms, this physics has been studied in a wide 
variety of systems ranging from driven nonlinear oscil¬ 
lators ET21 - fToj . superconducting circuits 133 ESI and cir¬ 
cuit quantum electrodynamics [.49] to phase transitions in 
driven-dissipative many-body systems [27] . It also sets a 
limit to the minimum mechanical temperature achievable 
using cavity cooling [3D]. 

Here, we show that quantum heating effects in an op¬ 
tomechanical system lead to observable signatures in the 
output spectrum of the cavity, i.e. the spectrum of out¬ 
put light that would be measured using a photomulti¬ 
plier. Unlike an OMIT experiment, such an experiment 
probes both the DOS of the cavity (as modified by the 
optomechanical interactions), as well as the effective tem¬ 
perature of the optomechanical polaritons. This “quan¬ 
tum heating” is already present at the level of the lin¬ 
earized theory of the optomechanical cavity; here, we in¬ 
vestigate how it is modified by the nonlinear interactions. 
Our work is timely, given that the generation of real pho¬ 
tons from mechanical vacuum fluctuations was recently 
probed experimentally by Lecocq et al. 1311 . 


A. Main findings 


Our work has many technical aspects to it that will 
hopefully be an aid in further studies of nonlinear quan¬ 
tum optomechanics (e.g. optomechanical lattices 0U- 
146] 1. It also predicts several new physical phenomena, 
in particular: 


1. Polariton thermalization 


In the standard linearized theory of a driven op¬ 
tomechanical cavity, quantum heating leads to very 
different effective temperatures for the two polariton 
modes (i.e. normal modes of the linearized theory), c.f. 
Eqs. (24). The nonlinear interaction tends to dilute this 
effect, as it allows energy exchange between the polari¬ 
tons and favours their thermalization. We study this 
competition in detail for two representative cases: a laser 
drive at the red-mechanical sideband (see Fig. [8]), and a 
laser drive detuned further to the red where nonlinear ef¬ 
fects are more important (see Fig. 10). The effects of the 
nonlinear interaction on this quantum heating physics 
can be seen experimentally in the cavity output spec¬ 
trum by tuning the nonlinear interactions into and out 
of resonance, see Fig. [9] 


2. New instabilities 

For a red-detuned laser, we find that a leading-order 
treatment of the nonlinear interaction suggests new kind 
of parametric instability not present in the linearized 
theory, where one polariton mode acts as an incoherent 
pump mode for the other, see Sec. |V D| Including higher 
order terms stabilizes the system as expected, see Fig. 0 


3. Two-phonon cavity heating 


We show that nonlinear heating of the cavity (as man¬ 
ifest in the cavity spectrum) are greatly enhanced for 
red-detuned laser drives near the second mechanical side¬ 
band (a detuning A ss —2%); this is a simple con¬ 
sequence of the lower optomechanical polariton being 
mostly phononic in this regime, causing effects to be 
enhanced by the typically large phonon lifetime (see 
Fig. [12]). 


B. Organization of the paper 


The remainder of this paper is organized as follow. 
We begin in Sec. [XT] by introducing the basic driven op¬ 
tomechanical system studied in this work, and review¬ 
ing the linearized theory. We also pay careful attention 
to the coupling to dissipative baths, and to how quan¬ 
tum heating effects can arise even without nonlinear¬ 
ity. In Sec. m we introduce the basic aspects of the 
Keldysh technique as applied to optomechanics, focus¬ 
ing first on the linearized system. In Sec. |IV[ we de¬ 
velop a Keldysh perturbation theory to treat the nonlin¬ 
ear optomechanical interaction, and introduce our self- 
consistent approach. In Sec. [Vj we discuss how the non¬ 
linear interaction modifies the non-equilibrium physics of 
the system. We present a physical picture where interac¬ 
tion effects can be mapped onto a couplings to additional 
“self-generated” dissipative baths. We also investigate in 
this section the new kind of parametric heating which 
arises with a red-detuned laser. In Sec. m we discuss 
the observable consequences of our predictions on the 
cavity output spectrum. We present our conclusions in 
Sec. EH] 


II. SYSTEM, TREATMENT OF DISSIPATION 
AND LINEARIZED THEORY 


A. Hamiltonian in the polariton basis 

We consider a standard optomechanical system where 
the frequency of a driven cavity mode is modulated lin¬ 
early by the position of a mechanical resonator (Fig. [I]). 
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The Hamiltonian governing its dynamics is given by mm 

H =u>cc^a + uiM^b + ga)a(b + E) 

+ i{a iVL e~ lULt d) - H.c.) + H diss . ( 1 ) 

Here, a is the cavity mode, with frequency wc and b 
the mechanical mode with frequency u>m- The parame¬ 
ter g is the single-photon optomechanical coupling and 
ci i n is proportional to the amplitude of a classical drive 
at frequency ujl- Finally, H d i ss describes dissipation due 
to the coupling to bosonic environments (both for the 
mechanics and the cavity). Going into a rotating frame 
at the drive frequency and displacing the cavity field by 
its classical value, induced by the coherent drive (i.e. 
a(t ) — > e~ lulLt (a + d(t))), the Hamiltonian of Eq. (jTJ) 
can be expanded into a quadratic part, known as the 
linearized optomechanical Hamiltonian, and a nonlinear 
interaction term such that 


11 = //i. + 11 


NL 


' Hdisst (2a) 

H l = -A^d + uj M b f b + G(d + d^)(b + b f ), (2b) 
H Nh = g$d(tf +b). (2c) 

In the above, we have defined the laser detuning A = 
«l- u>c and the many-photon coupling constant G = ag; 
we take a, g > 0 without loss of generality 55] ■ 

Our general approach is to diagonalize the quadratic 
part of the Hamiltonian, Hl (thus treating it exactly), 
and then treat the terms due to nonlinear interaction 
(Hnl) as a perturbation. Unlike treatments based on 
polaron-transformed Hamiltonians mm, we do not re¬ 
quire the coherent driving of the cavity to be so small 
that it too can be treated perturbatively. We focus ex¬ 
clusively on a red-detuned laser (A < 0), and only re¬ 
quire the cavity drive to be weak enough so that G 2 < 
— iv m A/4 = G 2 lit . For drives stronger than this critical 
value (or for a blue-detuned laser), the linearized coher¬ 
ent Hamiltonian Hl is unstable, and corresponds to a 
detuned parametric amplifier driven beyond threshold. 
This critical value also coincides with the onset of the 
well-known static optomechanical instability mm- Op¬ 
erating near this instability has been investigated as an 
alternative promising way to enhance the nonlinear inter¬ 
action m- Note that this instability is also closely anal¬ 
ogous to the superradiant phase transition in the driven 
Dicke model studied in m 


Focusing on A < 0 and G 2 < —com A/4, Eq. (2b) can 
be diagonalized via a Bogoliubov transformation to yield: 


11 \ — ^ ) llrr cjj l‘n 

<7 = ± 


(3) 


Here, c± destroys an excitation in the eigenmode of Hl 
with energy E± > 0, given by 

E± = ^ (“m + A 2 ± ^( W ^-A2)2-16G 2 Aw m ) ' ■ 

(4) 


Note that H_ tends to zero as G approaches the critical 
value at the instability G cr it- As Hl does not conserve 
the total number of photons and phonons, the operators 
c_, c+ mix photon/phonon annihilation and creation op¬ 
erators: 


b = a.b-C- + ctb,+c + + ab,-c_ -)- 0 ^+ 6 +, (5a) 

d = cxd-C- + a d ,+c + + a d -cl_ + a d , + c\. (5b) 


The coefficients ab/ d ,± and ctb/ d .± are functions of A /com 
and G/u m ; th eir explic it form in the case A = is 
given in Eqs. © ~( |A4[ ) of Appendix [Aj As the excita¬ 
tions described by c_, c + have both phononic and pho¬ 
tonic components, we refer to them as polaritons in what 
follows. 

Having diagonalized the linearized optomechanical 
Hamiltonian, we now express the nonlinear Hamiltonian 
Hnl in the polariton basis. We obtain interactions that 
do not conserve the total number of polaritons, 


#NL = ^ (ffaa'o-"cic f a'€" + 9aa’ + H.C.) 

~\~ (A_c_ -\- A- |_c_|_ H- H.c.). (6) 


A / D 

Here, the cons tant s g nn , n ,i and A a are all proportional 
to g (see Eqs. of Appendix [Aj). Note that the 

linear terms in Hnl (oc A a ) arise from normal ordering 
Hnl in the polariton basis; physically, while the zero- 
polariton state is the vacuum of Hl, this is no longer 
true when we include the nonlinear interaction. 


B. Resonant polariton interactions 


As discussed in previous works jlBI [IT] (and later in 
f-4'Sj). we can enhance the effects of even a weak single¬ 
photon coupling g by tuning G and A such that the non¬ 
linear processes in Eq. that scatters a + polariton 

into two — polaritons (oc gE _|_) become resonant. This 

requires E + [ A, G] = 2H^[A, G]; for a given laser detun¬ 
ing in the range A £ [— 2ojm, -wjif/2], this can always be 
achieved by tuning G = G res [A], where 


G res [A] = 


yi7AH-4(A 4 +0 

(icy— a W m ) 


(7) 


Once G is tuned to achieve this resonance condition, 
one can show using standard perturbation theory that 
all non-resonant nonlinear processes are suppressed by 
a factor of k/(E + — H_) oc k/lom compared to the res¬ 
onant process, where n is the cavity damping rate (see 
m for more details). In addition, in this regime, the 
relative modification of the polariton energies and wave- 
functions due to the linear terms in Eq. © are strongly 
suppressed by a factor of ( g/ujM ) 2 - Thus, if we focus 
on parameters near this resonant regime, in the resolved 
sideband regime (k/lom -C 1), and for weak nonlinear 
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interaction (j/w« <C 1), we can both ignore the renor¬ 
malization of polariton energies and wavefunctions, and 
make a rotating wave approximation on Eq. keep¬ 
ing only the resonant interaction. The coherent system 
Hamiltonian in this regime thus reduces to 

H e s = ^2 Eaclca + g{c\c-C- + H.c.), (8) 

<7 —± 


where g = g^ _, is the effective nonlinear coupling (see 

Eq. (A8) in Appendix [A|. Shown in Fig. |2](a) is the de¬ 
pendence of G res as a function of laser detuning A, as 
well as the behaviour of g on A, when G is tuned to be 
G res - For the rest of this paper, we focus on the dynamics 
governed by the effective coherent Hamiltonian H e g. 


C. Coupling to dissipative reservoirs 


We now turn to a more careful consideration of the 
effects of dissipation on our system. As is standard in 
optomechanics, both the mechanics and the cavity are 
taken to be linearly coupled to independent, Markovian 
bosonic baths (i.e. baths with constant DOS over the 
frequency range of interest). This is analogous to the 
approach taken in input-output treatment of dissipation 
(see, e.g. JT),2D!)- As we are interested in possibly large 
many-photon optomechanical couplings G, care must still 
be taken, as the eigenstates of our coherent Hamiltonian 
are polaritons, not individual photons or phonons; similar 
issues have recently been addressed in strongly-coupled 
circuit QED systems m ■ Moreover, the strong driving 
of the cavity can also lead directly to “quantum heating” 
effects, which manifest themselves directly in the treat¬ 
ment of the cavity dissipation. We describe these effects 
more in what follows. 


1. Cavity dissipation in the presence of optomechanical 
coupling and driving 

Consider first the interaction between the cavity and 
its dissipative reservoir. In the original lab frame, the 
cavity-bath interaction will have the following generic 
form gg 




FIG. 2: (Color online) (a) G les , the value of the many-photon 
coupling which leads a resonant nonlinear interaction, as a 
function of laser detuning A (c.f. Eq. 0k the correspond¬ 
ing effective nonlinear coupling g (c.f. Eq. 0) is also plotted, 
(b) and (c) Coefficients relating the optomechanical polaritons 
(eigenstates of the linearized Hamiltonian) to original photon 
and phonon operators (c.f. Eqs. |5|) as a function of laser de¬ 
tuning A with G = Gres [A]. In (b), we plot the “normal” co¬ 
efficients which relate polariton destruction operators to pho¬ 
ton/phonon destruction operators, i.e. ctb, a = (0, Ojcblr |0, 0), 
with 10,0) being the polariton vacuum. In (c) we plot the 


“anomalous” coefficients, i.e. on, }CT = (0, 0|c CT ft|0, 0); these co¬ 
efficients are directly related to the existence of “quantum 
heating” effects in the linearized theory (c.f. Eq. (15 0- 


drive frequency via the unitary 


£.=5>/7/i> 


Hl nt = i 


\k Pc kk f ’) {a+a,] ■ c==p 

-iu; L t ( ka + ^/j/j ] 

(9) 

L V i J\ 


( 10 ) 


where fj is the anihiliation operator for cavity bath mode 
j (frequency ujj ), re the damping rate of the photons in¬ 
side the cavity and pc is the bath DOS. As we consider 
a Markovian bath, we take re and pc to be frequency 
independent. 

We next transform to an interaction picture at the 


Note that we transform both cavity and bath oper¬ 
ators. The result is that in our interaction picture, 
terms in Eq. ([9]) which conserve excitation number are 
time-independent, whereas the remaining excitation non¬ 
conserving terms are rapidly oscillation at a frequency 
±2wz, ~ ±.2u>c (wc ^ | A|). Since up is much larger than 
all other energy scales in the rotating frame (i.e. ujc 
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E a , | A | ,uiM,g, k)i and since the bath oscillators necessar¬ 
ily have positive energies u>j > 0 , these rapidly oscillating 
terms can never become resonant. We can thus safely 
make a rotating-wave approximation, and drop them. 

With this rotating-wave approximation, the cavity- 
bath Hamiltonian then becomes time independent. Re¬ 
writing the cavity annihilation operator as per Eq. (2b I, 


one obtains a standard-rotating wave system-bath inter¬ 
action 

h k = (wj - u L )fjfj - 53 ( n ) 


H[ nt = i 


< 12 > 


Note crucially that in this final rotating frame, the 
transformed bath frequencies ujj can be negative (they 
extend down to — ul )■ While this may seem innocuous, 
things become more interesting when we re-write this 
interaction in terms of polariton operators (c.f. Eqs. ©): 

H™ 1 = fj ( a d,*c a + d d ^ci) + h.c. (13) 

' “ 3,° 


The anomalous terms which create or destroy two ex¬ 
citations here should not be dropped; as the polariton 
energies E a <C ~ u>c, these terms can be resonant 
in our interaction picture, as bath modes having Cjj < 0 
can be involved. Physically, such processes involve the 
creation of both a polariton and a bath excitation, while 
at the same time a (classical) drive photon is absorbed. 
Such heating processes thus involve the interplay of the 
system drive and the bath zero-point fluctuations, and 
are at the heart of quantum activation I.32l - I35| . 


2. Quantum heating and effective temperature 


The anomalous, excitation non-conserving terms in 
Eq. (13) can lead to polariton heating even if the cav¬ 
ity dissipation is at zero temperature. We can natu¬ 
rally associate an effective temperature to this heating 
by computing Golden Rule transition rates (50| . Con¬ 


sidering first the lower-energy — polaritons. Hf* will 
cause transitions uphill in energy between an initial state 
having N — 1 polaritons and a final state having N po¬ 
laritons at a rate Tjv.jv- l oc iV|<5:^,— | 2 - Similarly, it will 
cause transitions downhill in energy from the N to N — 1 
polariton state at a rate T n-i,n oc N\ctd-\ 2 - If these 
transitions were due to a bath in true thermal equilib¬ 
rium at temperature T, detailed balance dictates that 
Ti/r N—i,N = exp(— E_/k B T). In our case, we can 
use the ratio of these rates to define the effective tem¬ 
perature T£ av of the cavity dissipation as seen by the — 
polaritons: 




rjy,jy-i _ Q d,- 

T N -1,N o?d- 


(14) 


If these transitions were the only dynamics of the — po¬ 
laritons, they would indeed cause them to reach a thermal 
state at temperature Tf. av , with a thermal occupancy: 


n B [E_,T™ v ] = 


1 


d,- 


,E-/k B T™ 


- 1 


(15) 


x d,— 


We thus see the two crucial ingredients needed to ob¬ 
tain a non-zero effective temperature in a bosonic system 
where the physical bath temperature is zero. We needed 
both a coherent drive (yielding effective negative energy 
bath modes in our interaction picture), and a coherent 
parametric-amplifier type interactions (i.e. coherent in¬ 
teractions which do not conserve particle number, and 
hence yield dd- 7 ^ 0). In a completely analogous fash¬ 
ion, one can also associate an effective temperature T? av 
describing the quantum heating of the + polaritons. We 
stress that these effective temperatures have nothing to 
do with nonlinear interactions. 

For a concrete example of this effective temperature 
physics, consider the special case of a cavity drive at 
the red mechanical sideband, A = —u>m- One finds 


(c.f. Eqs. (A1)-(A4|): 


n B [E ± ,T : T] 


A=—ujm 


(l-y/l±2 G/u>m)‘ 
4-y/l ± 2G/um~ 


(16) 


Note that n B [E_,Tff v ] diverges when G approaches the 
onset of parametric instability, G —> ujm/ 2. In this limit, 
E_ —► 0, and the divergence of n B [E _, X^ av ] is equivalent 
to a fixed effective temperature T£ av ~ ujm/ 4 = |A|/4 
near the instability. Such behaviour is generic for quan¬ 
tum heating near parametric instabilities. 

The type of quantum heating phenomena described 
here is generic, and is found in a variety of related sys¬ 
tems, though the generic nature of the mechanism is of¬ 
ten not appreciated. Analogous effective temperatures 
also arise in mean-field treatments of driven dissipative 
phase transitions. For example, the study of the driven- 
dissipative Dicke model in Ref. [27] finds an effective tem¬ 
perature near the transition identical to that quoted af¬ 
ter Eq. (16). The physical origin is analogous to that 


in our system: it arises from the interplay of coherent 
parametric-amplifier interactions combined with a coher¬ 
ent linear driving. 


3. Mechanical dissipation 

We now turn to the interaction between the mechanical 
resonator and its dissipative bath. The starting interac¬ 
tion is analogous to Eq. © for the cavity dissipation, 

HT=iJ^Y{p)-9?j (b + tf), 

3 v j 

(17) 
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where gj is a bath annihilation operator, 7 the mechan¬ 
ical damping rate and pM is the constant DOS of the 
(Markovian) bath. In what follows, we consider the me¬ 
chanical bath to be at temperature Tm and define the 
mean number of excitations inside the bath at as 


n th = n B[WM, Tm], 


(18) 


with n B [co,T] being the Bose-Einstein distribution. 

As the mechanics is undriven, the bath Hamiltonian 
and bath-systenr interaction are unchanged under the 
transformation of Eq. (10 1 to the rotating frame at the 


laser frequency. This time, we re-write the interaction in 
terms of polariton operators using Eqs. ©> before making 
further approximations. This difference from the treat¬ 
ment of cavity dissipation stems from the fact that unlike 
u>c, wm is comparable to E a . It is thus crucial to go to 
the eigenstates basis of polaritons before assessing which 
terms may be safely dropped. In the polariton basis, we 
have: 


H\ nt =*, 


X ( 9 ] - 9 j) (a b ,a + a b ^) (c a + cl) . 


(19) 


We can now consider the role of terms in Eq. (19) that 
do not conserve excitation number. In contrast to our 
treatment of the cavity dissipation, here such anomalous 
terms can be dropped in a rotating-wave approximation. 
As there is no mechanical drive, there are no effective 
negative energy mechanical bath modes, and hence these 
terms can never be made resonant. Thus, we finally ob¬ 
tain a simple rotating-wave interaction between the me¬ 
chanical bath and the polaritons: 


m nt = 1 


7 


2np M 


X [(“" 


+ a a )g]c r7 + H.c. (20) 


j) 0 ’=± 


As there are no excitation non-conserving terms in 


Eq. (20), it is easy to confirm that the effective tempera¬ 


ture of the mechanical dissipation seen by the polaritons 
is simply equal to the physical temperature of the me¬ 
chanical dissipation. As already emphasized, quantum 
heating requires both the presence of coherent parametric- 
amplifier type interactions and the presence of a coherent 
drive. Here, while the coherent pararnp interactions are 
present, there is no coherent driving of the mechanics; as 
such, there is no quantum heating effects involving the 
mechanical bath. We see that even if both the cavity and 
mechanical baths have identical physical temperatures, 
the polaritons see them as having different effective tem¬ 
peratures. The driven nature of the system thus gives us 
interesting non-equilibrium physics even at the level of 
the linearized (i.e. quadratic Hamiltonian) theory. 

We end this subsection with a caveat on the valid¬ 
ity of treating dissipation via Markovian baths. For the 
cavity dissipation, this is an excellent approximation, as 
we are always probing the bath in a narrow interval of 
width ~ E a around up, an interval over which the bath 


DOS can be treated as constant (recall that E a <C wc). 
In contrast, a similar statement does not hold for the 
mechanical dissipation: we will be probing the mechan¬ 
ical bath at frequencies E a which could be significantly 
different from the mechanical frequency u>m- As such, 
it is not a priori obvious that the bath spectral density 
can be treated as flat. For simplicity we will nonetheless 
use the Markov bath approximation for the mechanics in 
what follows (consistent with the majority of works in op- 
tomechanics). For the weak dissipation limit of interest, 
the main effects of a non-flat bath spectral density could 
be easily incorporated into our calculations. One would 
simply make the mechanical contribution km to the in¬ 
trinsic polariton decay rates (c.f. Eq. (23)) proportional 
to the mechanical bath density of states at the relevant 
polariton energy. 


D. Lindblad master equation and effective 
polariton dissipation 

For further insight, it is useful to use the form of the 
system-bath couplings in Eqs. (12) and (20) to derive an 


approximate Linblad master equation for the dynamics 
of polaritons in our system. While we do not use such 
a master equation in our analysis, it provides a useful 
comparison point. Using the standard derivation (see, 
e.g., [49]) valid for weakly coupled Markovian baths, we 
obtain: 


p(t)=-i H e g, p(t) 


+ ^-(n® + 1 )D[c-]p + ^n°_D[cl_]p , 


+ ^(n° + + l)D[c + ]p+^n° + D[^]p, 


( 21 ) 


with the Lindblad super operator D[c a \p defined as 

D[c a \p = 2c a p(t)cl - clc a p{t ) - p(f)c£.c CT . (22) 

The effective polariton damping rates appearing in 
Eq. (21) are 


= 7 (a bjlT + a M ) 2 + k (a 2 da - a 2 da ) = + , 


(23) 


Here, we have introduced re“ ech («£ av ) as the contribution 
to the damping rate of the a polariton coming from the 
interaction with the mechanical resonator (cavity) dissi¬ 
pative bath. The corresponding effective bath thermal 
occupancies are 


h° = 


1 r 

H(J L 

^mech 


7 + a b}C ,) 2 n B [E a ,T M ] + na 2 d o 

n B [ E a , T m ] + K™n B [E a , T™ v ] 

i^mech 1 i^cav 


(24a) 

(24b) 


Here, Tm is the (physical) temperature of the mechanical 
bath (cf. Eq. ©), T^ av is the temperature of the cavity 
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bath as seen by the a polariton (cf. Eq. ©) and the 
different coefficients a are given in Eq. ([5]). As expected, 
Eq. (24b I represents a bosonic mode coupled indepen¬ 
dently to two disipative baths. An analogous expression 
holds for the effective mechanical occupancy used to de¬ 
scribe cavity-cooling experiments [40]. Note that we have 
made a standard secular approximation, allowing us to 
drop dissipative terms that do not conserve the number 
of each polaritons independently in Eq. (21 1 ; this is valid 
for the regime of interest E a , \ E + — E_ | n, 7 . 


The thermal occupation number of the effective baths 
given in Eqs. [24] and their corresponding temperatures, 
defined as n^E^T®] = n°, are plotted in Fig. [ 3 ] as a 
function of the detuning A. Anticipating our interest in 
nonlinear interactions, for each A we adjust the control 
laser amplitude so that G = G res (i.e. the value that will 
make the nonlinear interaction resonant, c.f. Eq. 0 ); 
this can be done for any A £ [— 2wm, - Wm/2]- One 
sees that even when the physical bath temperature is 
zero, quantum heating effects can yield effective polari¬ 
ton temperatures as large as ~ 0.1 quanta (solid curves). 
At non-zero physical temperature (dashed curves), these 
quantum heating effects persist, but are swamped by the 
contribution of mechanical noise at the edges of the de¬ 
tuning range considered. This is simply because near the 
limits of the detuning range, one polariton species is al¬ 
most all phononic (see Fig. [ 2 ]) and becomes very sensitive 
to thermal fluctuations of the mechanical bath. 


Returning to zero physical bath temperature, another 
striking feature in Fig. [3](c) is the sudden drop in effective 
temperature for the more phonon-like polariton branch at 
the edges of the detuning interval. As it will be of interest 
in what follows, we discuss this behaviour for detunings 
near A = —2 um in more detail; here, one sees a sudden 
drop in h°_ and T°. As A — >• —2 u>m, G res —>• 0, and the 
— polariton becomes simply a phonon. Expanding n°_ to 
lowest order in G res /iOM for Tm = 0 and 7 -C k, one gets 


n°_ = 


~ a d,~ 


1 g£ 
9 T? 


-K 


Of 4- ? 


8 G? e , 


(25) 


We see there are two competing effects associated with 
non-zero G res . The numerator reflects the parametric 
heating associated with the linearized optomechanical 
interaction. The denominator in contrast reflects that 
the phonon-like polariton has its lifetime decrease as 
G res increases and it becomes more photon like; this is 
just standard optomechanical optical damping. The re¬ 
sult is that the net quantum heating is maximized for 
[G re s/wM ] 2 ~ 7 /k), corresponding to a laser detuning 
A + 2 lom ~ 7 /re; the maximum occupancy h° that can 
be obtained is | (as seen in Fig. §b))- 



FIG. 3: (Color online) (a) Damping rates of the a = ± 
polaritons (cf. Eq. [23]) in the linearized theory (i.e. g = 0), as 
a function of the detuning A. For each A, we tune the cavity 
control drive to maintain G = G res (A). Near A = —2c Om, the 
- (+) polariton damping tends to 7 ( k ) as it is mostly phonon 
(photon) like; the converse is true near A = — wm/2. (b) 
Thermal occupation numbers h° of the effective bath coupled 

l) also when g = 0 and G = 


to the a polaritons (cf. Eq. (24 
G res . Here, characterises t 
the mechanical bath (cf. Eq. 


le (physical) temperature of 


18 1 ). (c) The temperatures 


corresponding to the thermal occupation (fee = 1 ) in the 
same regime as panels (a) and (b). For each curves, 7 /re = 
10 ~ 4 and uim/k = 50. 


III. KELDYSH DESCRIPTION OF THE 
LINEARIZED SYSTEM 


As shown in the previous section, the driven nature 
of the system leads to non-equilibrium physics even at 
the level of the linearized theory (i.e. the two polari¬ 
ton species see different effective temperatures). Con¬ 
sequently, we need to use the Keldysh formalism 2S] in 
order to describe the dynamics and to properly construct 
a perturbation theory that treats the nonlinear interac¬ 
tion present in H e g (Eq. ^). In this section, we quickly 
introduce this approach, considering first the linearized 
Hamiltonian of Eq. (2b) and the couplings with the en- 






















vironment given by Eqs. (12) and (|20[). We will also map 


the resulting free Keldysh theory onto the simple master 
equation previously derived (Eq. ©), again working in 
the regime where polaritons have a well defined energy 
(i.e. n a <C E a ). We stress however that even without 
interactions, the Keldysh approach is more general than 
a Lindblad-style master equation, as it is not restricted 
to Markovian baths. 

In the Keldysh formalism, we represent our linearized 
optomechanical system by a field theory which is general 
enough to allow the system to be in an arbitrary, non¬ 
equilibrium state. In this field theory, there are two time- 
dependent fields (classical and quantum) corresponding 
to each annihilation operator in the original theory. Con¬ 
sequently, the quadratic action that conserves the num¬ 
ber of particles of our two independent bosonic modes 
(polaritons) will have the following general form: 




—oo —oo 


(26) 


used in our previous work DU- Here, we instead follow 
a different but equivalent route to obtain the bare Green 
functions and Keldysh action. 

The goal here is to write the different Green functions 
such that our description of the linearized theory in the 
Keldysh formalism is completely equivalent to the mas¬ 
ter equation of Eq. pi]). As discussed before, Eq. pl|) 


describes two independent bosonic modes with simple 
Markovian damping rates n a (cf. Eq. (231) and energies 


E a (cf. Eq. 0). Consequently, the retarded (advanced) 
Green functions, which give the response functions of the 
system, have to adopt the simple following form: 


g£m = {g£m} = 


w - E a + iK'j/2 ’ 


(30) 


i.e. a polariton has an energy E a and a lifetime 1 /k ct . 

We now construct the Keldysh Green functions using 
the same approach. As is standard [23], we define a dis¬ 
tribution function f[uj] that relates the Keldysh and the 
retarded Green functions, such that 


G a M = -2i(2f[u] + l)Im[G~[w]]. 


(31) 


Here, the c aq / c i(t) are complex functions of time and 
G” 1 (t — t!) is the (operator) inverse of the unperturbed 
(i.e. g = 0) Green function. The latter is given by a 2 x 2 
matrix 


G a {t) 


G?(t) G*(i)\ 

0 




(27) 


In terms of Heisenberg picture operators, each element is 
defined as 


G«(t) = {Gait)}* ee -iOit)([cait),c t (0)]>, (28a) 

G*(t) = -*<{<■,(*),$ t(0)}>, (28b) 


where the expectations are taken with respect to the ini¬ 
tial density matrix without nonlinear interaction (<7 = 0 ). 
Here, G*(t) and G${t) are the standard unperturbed 
retarded and advanced Green functions, which govern 
the linear response properties of the unperturbed sys¬ 
tem. They are also related to the unperturbed DOS of 
each polariton, given by 


P°M = —Im [G«M: 
7r 


(29) 


Finally, it) is known as the (unperturbed) Keldysh 
Green function. It encodes knowledge of the energy dis¬ 
tribution function of each polariton (as we will see more 
clearly below). 

With these definitions in hand, one could follow the 
standard approach used in input-output theory and de¬ 
rive the Heisenberg-Langevin equations from the coher¬ 
ent Hamiltonian of Eq. ([ 8 ]) and the particular form of 
the system baths coupling given by Eqs. (12 1 and (20 1 . 
From there, one can directly get the bare Green func¬ 
tions by calculating Eqs. (28 1 . This approach has been 


This function f[ui] parameterizes the occupation of dif¬ 
ferent polariton energy eigenstates. As an example, for a 
system in thermal equilibrium at a certain temperature, 
f[uj\ would be the corresponding Bose-Einstein distribu¬ 
tion and Eq. (31) would be an exact statement of the 
fluctuation-dissipation theorem [25] . In the particular 
case studied here, the free polaritons have sharply peaked 
single particle DOS p a [ui\ (cf. Eqs. (29) and ( |30| with 
E a ), such that for the u polariton, the function 
f[uj] of Eq. (31) can be approximated as /[w] ~ f[E a ]. 
Finally, if we insist that the average occupancy of the po¬ 
lariton matches that in the master equation description, 
then we must have f[E a \ = n°. We thus have 


G^[ W ] = -2i(2n“+l)Im[G"M]. 


(32) 


Describing the linearized theory in the Keldysh for¬ 
malism using the bare Green functions (30) and (32) is 


thus completely equivalent to the Lindblad master equa¬ 
tion (21). We recall that the two assumptions underlying 


these two equivalent descriptions are that the polaritons 
have sharply peaked DOS (i.e. dissipation is weak) and 
that coupling between the + and — polaritons due to dis¬ 
sipation is negligible (i.e. secular approximation made to 
derive Eq. pi ), see Sec. |II Dh. 


Finally, we note that it is possible to derive exact 
Langevin equations from the linear Keldysh action given 
in Eq. (26). Briefly, one first decouples the quadratic 


quantum-field terms via an exact Hubbard-Stratonovich 
transormation; this introduces new fields £® nv (t) which 
have a Gaussian action. One can then exactly do the 
integrals over quantum and classical fields. The result¬ 
ing functional delta function corresponds to the Langevin 
equations: 

d t Cait) = - (lEa + ’f) Cait) - £" V (i), (33) 
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where the noise £® nv (f) is Gaussian with zero mean. The 
only non-zero noise correlation functions are given by 


(£ m (t )[£' v (OH = + 1/2 )S(t - t')<W- (34) 


As expected, Eqs. (|33j) and (341 represent two uncoupled 
damped harmonic oscillators each in contact with their 
respective finite temperature Markovian baths. 


IV. KELDYSH PERTURBATIVE TREATMENT 
OF POLARITON INTERACTIONS 



A. Self energies and dressed Green functions 



Have established the Keldysh formulation of the lin¬ 
earized optomechanical theory, we can now address the 
effects of the nonlinear interaction as a perturbation. We 
assume throughout this section that the drive laser has 
been tuned to make the non-linear interaction resonant, 
i.e. G = G res , and thus consider only the resonant inter¬ 
action process given in Eq. Q. 

First, the action generated by the nonlinear interaction 
in H e s (Eq. in the cl-q basis is 


5*NL — r— j dt ^C_|_ gC_ )C £C_ jC £ T 2c_j_ c iC— ^qC— 


V2J- 


Tc_j_ qC— ^qC— _q ~\~ C.C.) , 


(35) 


where the time dependence of the fields are implicit for 
clarity. Diagrammatically, each terms in the nonlinear 
action corresponds to a vertex shown in Fig[4|a). The 
vertices with a single quantum field correct the classical 
saddle point, 


d{S\ J + 5*nl) 

dc *+,q(t) 


C =t ,q —0 


9(5l + S NL ) 


8c*_Jt) 


C±,,=0 


r°° i 

/ dt’ [G+(t - 1')] c +jC ;(f') 

J — OO 

+ -^c_ ci (f)c_, c ;(f), (36a) 

/ OO 

dt’ [G^(t - t')] 1 c_ iC j(t') 

-OO 

+ V2gc + ^ c i{t)c*_ cl (t). (36b) 


and thus correspond to a classical nonlinear potential. 
The term with three quantum fields is more of a purely 
quantum effect. It could be interpreted as an effective 
nonlinearity of the quantum noise. 

One can still derive Langevin equations from the re¬ 
sulting nonlinear action if one ignores the terms which 
are cubic in quantum fields. In this approximation, one 
obtains modified versions of the Langevin equations in 
Eq. (|33|): 


dtC- = —(iE_ + —)c_ — 2 igc* + C- — £® nv , (37a) 
d t c+ = -( iE+ + / y)c+ - igc-c- - £+ nv , (37b) 


FIG. 4: (Color online) Diagrams representing the self energies 
up to second order in g in the Keldysh formalism. The clas¬ 
sical vertices are composed of only one quantum component 
while the quantum vertex is composed of three. 


where the autocorrelation functions of the ( t ) noise are 
unchanged by the nonlinear interactions. We have sup¬ 
pressed the explicit time-dependence of fields here for 
clarity. 

We will not use these approximate quantum Langevin 
equations further, but proceed in a way that does not 
neglect terms that are cubic in the quantum fields. As 
we are interested in weak nonlinear couplings g, we will 
compute the self energy (E[w]) of our Keldysh Green 
functions perturbatively to order g 2 . At this order, all 
relevant scattering processes (see Fig. [4]) conserve the 
number of polaritons independently (i.e. the self energies 
are diagonal in the +/- index). Consequently, the Dyson 
equation that gives the Green functions in presence of 
interactions can be separately written for each polariton. 




— GAud 1 


M £>]\ 
WM 0 J ' 


(38) 

Here, we used Q^' R ’ K [oj\ to distinguish the full Green 
functions (i.e. including the effects of g) from the unper¬ 
turbed ones G^ ,r,k [oj\. 


1. Retarded self energies and interaction-induced polariton 
damping 


The diagrams related to the second order retarded self 
energies are shown in Fig. |4jb) . From these, one straight¬ 
forwardly calculates 


E -M =c e -^f 

=cf^± 


+ K+ 

{E+ -E_) + i^=±^’ 

K- 

uj — 2 E_ + in_ ’ 


(39a) 

(39b) 


with E^[w] = {E^[u;]}*. As discussed in detail in 17j . 
the self energies describe the hybridization be¬ 

tween the near-resonant |+) and | —) polariton states. 
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The resonant nonlinear processes underlying this hy¬ 
bridization (see Fig.[4|b)) are responsible for the sharply 
peaked self energies of Eqs. (391 (k ct <C E a ). We have 


introduced effective cooperativites to parameterize 
how strong the decay rates resulting from these processes 
are (i.e. imaginary part of the self energy) on resonance, 
compared to the intrinsic polariton linewidth. Defining 


C t M = -2Im[SfM] 


we have 


ceff = r^[E + - E _] _ 16g 2 (n°_ - n° + ) 


C_(k_ + K + ) 


„eff _ ^[2E_] _ A~g 2 (2n°_ + 1) 
c + _ — 


K+ 


K-K+ 


(40) 

(41a) 

(41b) 


The definitions of Cf. ff ,C5® are analogous to the defi¬ 
nition of the standard optomechanical cooperativity C = 
4 G 2 /k7 as the cavity-induced “optical damping” of the 
mechanics to the intrinsic mechanical damping. The ef¬ 
fective cooperativities are plotted in Fig. [5] as a function 
of the detuning (keeping G = G res for all detunings). 

A crucial feature of the interaction-induced polariton 


damping described by Eqs. (41) is their explicit temper¬ 
ature dependence. This is a direct consequence of the 
multi-particle nature of the relevant decay process. For 
the — polariton, we have T 1 "* oc (n°_ — n+), as expected 
for a bosonic polarization bubble; a similar damping rate 
is found for an oscillator coupled quadratically to an os¬ 
cillator bath [52] . In true thermal equilibrium, the fact 
that E + > £L guarantees this factor is positive, yield¬ 
ing r 1 "* > 0. Our system however is not in true thermal 
equilibrium: as discussed in the previous section, it is 
possible to have (n°_ — n\) < 0 by having suitably differ¬ 
ent effective temperatures for the two polariton species. 
One thus finds that the interactions can lead to nega¬ 
tive damping: r™* < 0. The physics of this regime and 
the possibility of true instability are discussed further in 
Sec. [TUI 

The retarded polariton self energies presented here di¬ 
rectly lead to an interaction-induced modification of the 
polariton DOS, p a ]d\, given by 


p a [u}] = ~-Im [G„[<J\] = ——Im 

7T 7T 


1 


U-E a + in a /2 - [w] 

(42) 

Signatures of g in the DOS (and corresponding changes 
to OMIT-style experiments) were the focus of our pre¬ 
vious work m- For strong enough nonlinear coupling 
g , the self energies turn from simply describing an extra 
broadening of the polaritons, to describing the coherent 
hybridization of the resonant |+) and | —,—) polariton 
states. The key consequence of this is that the single peak 
in the + polariton DOS splits. Further details about this 
splitting (and how it can be measured via an OMIT-type 
experiment) can be found in Ref. 1171 


2. Keldysh self energies 

We now turn to the Keldysh self energies, which are 
also directly calculated from the diagrams of Fig. |4|c). 
Here, we parametrize each Keldysh self energies via a 
thermal occupancy factor n™* associated with the inter¬ 
action, defined such that 


EfM = -2i(2< nt + l)Im [E«M]. 


(43) 


This parametrization is always possible if we let the n™* 
to be frequency dependent. However, we find them to be 
frequency independent and given by 


(n° + 1) 


n_ = 


(n°_) 2 
2n°_ +1' 


(44) 


We stress that these results (as well as the self-energy 
results above) are based on only keeping the nonlinear 
polariton interaction in Eq. ([8]), and thus assume that 
E + ~2£L. 


B. Self-consistent calculation 

The self-energy results discussed so far (and in Ref. 1171) 
only retain diagrams to leading order in g. To capture 
higher order effects and effectively resum diagrams at 
all orders in perturbation theory, one can make the dia¬ 
grams in Fig. [4] self-consistent. One simply replaces all 
internal propagators in the diagrams by full dressed prop¬ 
agators. The self energy thus becomes a functional of the 
full dressed Green function, and the Dyson equation be¬ 
comes a self-consistent equation for the full Green func¬ 
tion. Solving this self-consistent Dyson equation allows 
to capture a particular ensemble of processes at all or¬ 
ders in g. Note that a related self-consistent Keldysh 
approach was previously used to study a nonlinear para¬ 
metric amplifier near threshold [2H - I2H) . 

In practice, one solves such self-consistent equations 
iteratively: in each step, one calculates the self energy 
using the current versions of the full Green functions, 
and then uses these to update the Green functions which 
will be used for the self-energy calculation in the next 
iteration. Applying this iterative procedure until con¬ 
vergence solves the self-consistent Dyson equation. To 
improve the accuracy of our approach, we have applied 
this iterative strategy for the results shown in Fig. [5] to 
Fig. |T2} For each calculation, we have performed 20 it¬ 
erations, which turns out to be more than enough to get 
excellence convergence of the Green functions. Note that 
by using this self-consistent approach, n™ 1 , as defined in 
Eq. (43), becomes frequency dependent. 


As expected, the self-consistent approach does not con¬ 
verge if the nonlinear interactions become too strong. 
Convergence is not solely controlled by the magnitude 
of g , but is primarily determi ned by the effective cooper¬ 
ativities C° s defined in Eqs. (41). These cooperativities 
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FIG. 5: (Color online) Solid lines: Leading-order-in-g effective 
cooperativities C^ s associated with the nonlinear interaction 
(defined in Eq. (411), as a function of the detuning in the 
resonant regime (G = G res [A]). One sees that the are 
enhanced near A = —2 lom- The inset zooms on the region 
where goes below zero, which signals the possibility of 
a new kind of instability, as discussed explicitly in Sec. |V D| 
Circles: results of the self-consistent perturbation theory de¬ 
scribed in Sec. |IV B| which includes diagrams at all orders in 
g. For all values of C)) ff shown here, the approach converges 
and is in good agreement with n ume rical simulations of the 
Lindblad master equation (cf. Eq. (211). The parameters used 
are y/K = 10~ 4 , g = k, u>m = 50k and = 0. 


involve both the magnitude of g and the polariton oc¬ 
cupancies, reflecting the fact that large temperature can 
also enhance the importance of nonlinearity. As shown 
in the figures, we find that when the self-consistent ap¬ 
proach converges, it also is in excellent agreement with 
full numerical simulations of the Linblad master equation 
describing the system, Eq. (21). 


V. INFLUENCE OF POLARITON 
INTERACTIONS ON NON-EQUILIBRIUM 
EFFECTS 


A. Interaction-induced effective environment 


From our previous discussion, we see that the non¬ 
linear interaction gives rise to self energies which mod¬ 
ify the single-particle properties of polaritons. The 
imaginary part of the retarded self energies describe 
new interaction-induced polariton damping rates 
(c.f. Eq. ©), whereas the Keldysh self energies describe 
new interaction-induced heating effects, with associated 
thermal occupancy factors h™ 4 (c.f. Eq. (44)). This sug¬ 
gests that in terms of single-particle properties, the effects 
of interactions are equivalent to having coupled the lin¬ 
earized optomechanical system to new dissipative baths. 
In what follows, we make this picture of an “interaction- 
induced effective environment” explicit. 

First, note that all single particle polariton properties 


of our system are described by the effective quadratic 
action 


Sf = 


oo 

JJ dtdt'[c* a c i{t), (f (r q {t))Q~ 1 {t — t') 


C-(7,cl{t' ) 
Ca.qit 1 ) 


(45) 


where Q a 1 (t — t') is the Fourier transform of the inverse 
of the 2x2 matrix of the full Green functions given in 
Eq. ((381). 


Further, as discussed in Sec m this quadratic action 
is completely equivalent to a set of linear Langevin equa¬ 
tions. Using the standard derivation [23], the effective 
action in Eq. (45) is equivalent to the Langevin equa¬ 
tions 

d t c a {t) = - (iE a + -^-)c a (t) - / dt'Y>K{t-t')c a {t') 

" J —OO 

(46) 

Here, E^ (t — t') is the Fourier transform of the retarded 
self energies given in Eqs. (391. The noise functions 
£® nv (t) and £” lt (f) describe complex independent Gaus¬ 
sian noise processes with zero mean. £® nv (t) describes the 


intrinsic polariton dissipation as discussed in section |HI[ 
its correlators are given in Eq. ( |34| ). The only non-zero 
correlator of the new noise (t)~is 

<e nt wk? con = i t‘(* - o« nt + i/2)<w. (47) 

where r™ t (t — t') is the Fourier transform of 
the frequency-dependent interaction-induced polariton 
damping given in Eq. (40). 


These Langevin equations reproduce the intuitive pic¬ 
ture sketched above: each polariton species is now ef¬ 
fectively coupled to two independent dissipative environ¬ 
ments, with corresponding damping rates K a and rj) 1 * [w], 
and corresponding thermal occupancies nff and hf [t . The 
first bath corresponds to intrinsic dissipation (i.e. the in¬ 
trinsic mechanical and cavity dissipation), whereas the 
second is due to polariton-polariton interactions. 

It is worth stressing that the induced damping rates 
r™ t [w] are in general sharply peaked functions of fre¬ 
quency, due to the resonant nature of the relevant scat¬ 
tering process. For some parameters, the width of F^ lt [cu] 
can even be much smaller than the width of the density 
of state p a [u\] e.g. the + polaritons for A near —2 ujm- 
In contrast, there are other cases where the density of 
state is much sharper than the interaction-induced dissi¬ 
pation rate, as is the case for the + polariton for A near 
—u>m/ 2. As a result, the “interaction-induced” baths 
cannot always be considered as Markovian. 

The thermal occupancies h’f 1 associated with the 
interaction-induced environments are plotted in Fig. [6] 
as a function of the detuning A, in the interesting case 
where all intrinsic dissipation (i.e. mechanical bath, cav¬ 
ity bath) are at zero temperature. 
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FIG. 6: (Color online) Effective thermal occupancies for the 
various baths coupled to the two polariton species, as a func¬ 
tion of A with G = Gres [A]. Solid lines: occupancies associ¬ 
ated with t he “ interaction-induced” baths to leading order in 
g, c.f. Eqs. (441. Circles: same, but calculated using the all¬ 
orders self-consistent approach of Sec. |IV B| Dashed lines: oc¬ 
cupancies of the intrinsic polariton baths, c.f. Eqs. ( |24[ ). Note 
the leading-order occupancy for the — polariton interaction 
bath (nT 4 ) diverges when n°_ = n+ (A/uom ^ —0.8) and be¬ 
comes negative when > h c i (A /ujm > —0.8), as shown in 
the inset. This divergence persists in the self-consistent the¬ 
ory, but occurs at smaller-magnitude detunings. This neg¬ 
ative occupancy signals the possibility of a new kind of in¬ 
stability, as discussed explicitly in Sec. |V D| All curves are 
plotted for = 0 (cf. Eq. (181), ■y/n = 10 -4 , ojm/k = 50 
and g = k. 


B. Interaction-induced quantum heating 


We now discuss in more detail the behaviour of Eq. (44) 
which gives the thermal occupancies h™ 4 of the effec¬ 
tive “interaction-induced” dissipative baths introduced 
in the previous subsection. For simplicity, we focus on 
the case of exact resonance, where G = G res and hence 
E + = 2 E_. Consider first the case where the linear- 
theory polariton dissipation is in thermal equilibrium at 
temperature T eq , i.e. = ns[£ ff ,T eq ]. In this case, it is 
easy to confirm that for each polariton, n“ 4 = n°, i.e. the 
“interaction-induced” dissipation also corresponds to the 
same temperature T eq . Thus, if without interactions the 
polaritons start in equilibrium at the same temperature, 
then the same is true with interactions. 

The actual situation is however more complicated: due 
to quantum heating effects, the effective temperatures of 
the two polariton species are different even without in¬ 
teractions. h°_ and fi° are thus not related as they would 
be in thermal equilibrium; this can be parameterized as 


n°, = 


ff!=G + 8n\. 

2 n°_ + 1 + 


(48) 


Thermal equilibrium and the condition E + = 2EL would 
imply Sn [{_ = 0; 6nff 7^ 0 means that even in the lin¬ 
earized theory, the two polaritons experience different ef¬ 
fective temperatures (see Fig. [3]). 


Using this definition, the thermal occupancy of the — 
polariton interaction-induced bath becomes: 


n” 4 =n°_+ Sn 




(1 + 2n°) 2 


(1 + n°_) — <5?r)[_(l + 2 n°_) 


(49) 


Thus, a deviation from true thermal equilibrium in the 
linear theory (i.e. without polariton interactions) causes 
the occupancy of the interaction-induced bath h™ 4 and 
the intrinsic bath h° a (linear-theory dissipation) to de¬ 
viate from one another. This is not surprising: in this 
case, the nonlinear interaction between the two polariton 
species tends to favour their thermalization, and hence 
transfers energy from the high-temperature species to the 
low-temperature species. 

Finally, we also stress that even in the case where 
the intrinsic mechanical and cavity dissipation is at zero 
temperature (i.e. the system only experiences vacuum 
noise), the interaction-bath thermal occupancies h™ 4 will 
be non-zero, and are in general different from n”.. This is 
shown explicitly in Fig.[6j We thus see that interactions 
change the effective temperature associated with quantum 
heating effects. 


C. Polariton energy distribution functions 


The picture established so far is that our optomechan¬ 
ical polaritons are each effectively coupled to two in¬ 
dependent effective environments, one of which is self¬ 
generated and due to the nonlinear optomechanical in¬ 
teraction. In the limit where both the intrinsic cavity 
and mechanical dissipative baths are at zero tempera¬ 
ture (n^ = 0), both these effective environments describe 
quantum heating physics. Together, they will determine 
the total number of polaritons produced by quantum 
heating, and more specifically, the energy distribution 
function of the polaritons. We define this distribution 
function in the standard manner, as an energy dependent 
distribution function [w]. This quantity is defined via 
the full Keldysh and retarded polariton Green functions: 

G ?M = -2i(2nf[w] + l)Im [Q?[w\] (50) 


If our polaritons were in thermal equilibrium at temper¬ 
ature T eq , then the distribution function would 

simply be the Bose-Einstein distribution ns[w,T eq ]. In 
contrast, in our system this function will be determined 
by the thermal occupancies of the two effective baths, 
and the strength of the couplings (i.e. damping rates) to 
each. Using the expression of the dressed Green functions 
Q[ui\ coming from the Dyson equation (cf. Eq. (38)) and 
the relation between the self energies given in Eq. (431, 
one finds 


„ M*. 


+ n a 


(51) 


This is exactly the simple expression that would be ex¬ 
pected for a free bosonic mode coupled independently 
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FIG. 7: (Color online) (a) Solid lines: net polariton occu¬ 
pancies in the presence of nonlinear interaction, nf[E a ], cal¬ 
culated using leading-order self-energies, as a function of A 
with G = Gres [A]. Circles: same, but calculated using the 
self-consistent approach (cf. Sec. IVB|. Dashed curves: Oc¬ 
cupancies for <7 = 0, i.e. calculated in the linearized theory. In 
each case, the mechanical temperature is zero (nj£ = 0 ) and 
g = k. (b) Same as (a), but now with = 100 and g = 0.1k; 
main plot and inset show different ranges of A. Near —2 lom, 
one sees an important contribution to the mean number of 
+ polaritons due to nonlinear interaction; this contribution is 
greatly enhanced by temperatures, as discussed in Sec. | VI D| 
All the curves are plotted for 'y/n = 10~ 4 and c um/k = 50. 


to two baths; the same form holds for the linear theory 
(cf. Eqs. @). 

In order to focus our attention on the contributio n of 
the nonlinear interaction to n® ff [w], we re writ e Eq. fl5l| ) 
using the expressions for T^w] (cf. Eqs. (39) and (40)) 


and nf (cf. Eq. 44). Doing so, one gets 


E M =n° +I a 


(u> — l 


+ 7i 


(52) 


with 


resonant nature of the relevant nonlinear scattering pro¬ 
cess. 

For exact resonance (E + = 2 E_), both n® ff [w] and 
the single particle DOS p a [ cu] (cf. Eq. (42)) are peaked 
at E a , so that the nonlinear interaction heating effects 
are maximal. Even in this case though, the frequency 
dependence of the interaction contribution to nf [w] can 
be very different than that of the polariton DOS. In this 
fully resonant case, the polaritons distribution functions 
evaluated at ui = E a adopt the following simple form 


Cfnf + 


^effrp 1 _ '-'<7 n c 1- 'V 
l£y - «■ +1 


(54) 


Eq. (54) is plotted in Fig. [?] as a function of the laser 
detuning A and is compared to nf \E a ] o btained using 
the self-consistent approach (cf. Sec. IVB). 


D. Nonlinear parametric heating 


Among the more striking non-equilibrium behaviours 
possible in the linear theory is the possibility of having 
n+ > n_, i.e. the thermal occupancy of the higher energy 
+ polariton exceeds that of the — polariton. We discuss 
this regime in more detail here, focusing on the exactly 
resonant case where E + = 2 E-. 

We start by recalling that the total damping of the — 
polariton is 


*_ ot [£;-] = k- + ri"*[£_] = «_ (i + cf). 


(55) 


From Eqs. (41), we see that if we have the occupancy 
inversion n[ J_ > n °_, then Cf < 0, and hence the contri¬ 
bution of the nonlinear interaction to the damping rate 
of the — polaritons becomes negative. This is at first 
glance surprising: we have opened a new scattering pro¬ 
cess for the — polariton via the nonlinear interaction, and 
yet we get an increase in its lifetime. We also have the 
possibility of an instability if Cf < — 1. 

From Fig. [5j one sees that a negative Cf occurs for 
detunings near A = —0.5 wm- In that regime, the — po¬ 
lariton is mostly photonic and the + polariton is mostly 
phononic. Consequently, by having a high intrinsic me¬ 
chanical bath temperatures nj^ 1, one naturally can 
achieve the inverted occupancy regime where n° + n°_. 

In that case and in the limit where 7 <C k, one can sim¬ 
plify the instability condition, i.e. Cf < — 1, to 


r>eft 

4 = (nf - n°) - + a CeS , (53a) 

7 - = K ~^ K+ yJl + Cf, 'r+ = K- y /l+Cf, (53b) 
w_ = E + — E-, = 2 E-. (53c) 

The contribution from the interaction-induced environ¬ 
ment appears as a sharp Lorentzian in the polariton dis¬ 
tribution functions. This is a direct consequence of the 


> n°_ + 


V(ft + K+) 


16 g 2 16g 2 ’ 

and the final mean number of — polaritons becomes 


(56) 


nf[E_] 


16<) 2 n+ 


k z _ — 16 <? 2 n+ 


(57) 


The surprising negative damping occurring here can 
be understood as the result of a parametric instability 
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arising directly from the polariton interaction. In fact, 
Eqs.(|56|) and (57) have exactly the same form as the in¬ 


stability condition and the mean number of excitations 
that one would find for a degenerate parametric ampli¬ 
fier (DPA) pumped near degeneracy [M1IM]. In a DPA, 
a pump-mode photon scatters into to signal mode pho¬ 
tons, and the pump mode is coherently driven. In our 
system, the — polariton plays the role of the signal mode 
in a DPA, while the + polariton plays the role of a pump 
mode that has been incoherently driven by noise. Despite 
this incoherent driving, the form of the above equations 
is the same as a coherently-driven DPA (see Appendix 
©• Note that non-degenerate parametric amplifier in¬ 
stability can be realized in a linearized optomechanical 
system driven with a blue-detuned laser, see e.g. j53U54]. 
In contrast, the instability described here occurs for a 
red-detuned drive. 

Note that our discussion here is based solely on using 
the leading order results for the polariton self energies. 
Including higher-order effects via our self-consistent ap¬ 
proach can dramatically change the onset and magnitude 
of the interaction-induced negative damping. We discuss 
this more in Sec. IVI Cl 


VI. OBSERVABLE SIGNATURES OF 
QUANTUM HEATING EFFECTS 

In the previous sections, we have demonstrated how 
quantum heating effects can lead to a finite density of 
optomechanical polaritons at zero temperature; we also 
discussed how these quantum heating effects can be mod¬ 
ified by the nonlinear interaction. In this section, we dis¬ 
cuss how these effects lead to observable signatures in the 
light leaving the optomechanical cavity. We first relate 
the cavity output spectrum to the polariton distribution 
functions, and then discuss specific parameter regimes 
where the heating effects are most prevalent. In addi¬ 
tion, we propose a way to effectively control the strength 
of the nonlinear interaction in experiments by tuning in 
and out the resonance condition (i.e. by varying G at 
fixed detuning A). Doing so, one can explicitly isolate 
and observe the nonlinear interaction signatures in the 
cavity output spectrum. 


A. Polariton energy distribution functions 

To measure polariton occupancies, we consider a mea¬ 
surement of the flux of photons leaving our cavity (assum¬ 
ing a single sided cavity, and that the reflected classical 
drive tone is filtered away). The spectrum of this flux is 
given in the standard manner [49] by the normal-ordered 
cavity spectrum (also known as the “lesser” Green func¬ 
tion within the Keldysh technique), 

/•OO 

S d [u] = / dte iut {S(0)d(t)). (58) 


Re-writing Eq. (58) in terms of polariton Green func¬ 
tions yields 


SdM =2 tt Y a| CT < ff Mp CT [u/| 

<T=± 

+ 2n Y v-lAnfi-u] + 1 )At[-w]- (59) 

CT —± 


Here, the coefficients and did,a are the change- 

of-basis coefficients introduced in Eq. (5a) and plotted 
in Fig. 2j and the polariton DOS, p a [(jj ]7 is defined in 
Eq. (|42). We are still working in a rotating frame with 
respect to the laser drive frequency wl, hence w = 0 im¬ 
plies output photons leaving at the laser frequency. The 
negative frequency term means that removing a photon 
at frequency w can involve creating a polariton at fre¬ 
quency —w, as expected from the presence of “anoma¬ 
lous” terms in Eq. (5b). Finally, Eq. (59) reflects the 
fact that the only non-zero polariton Green functions are 
those that conserve the number of polaritons indepen¬ 
dently, as discussed in Secs. m and |IV| 

Using the same assumptions, we also derive the cavity 
DOS and its energy distribution function: 


PdM = Y 

<T —± 


1 S d [v 

2n Pd [w] 


a d,*P<?[-u]) , 


(60) 

(61) 


From the cavity energy distribution, we can use the Bose- 
Einstein distribution to define an effective cavity temper¬ 
ature (ks = 1), 


Tf M 


UJ 


In 




(62) 


which is always possible if we let the effective tempera¬ 
ture to be frequency dependent. 

As discussed extensively in ca, the polariton DOS 
p CT [w] can be directly measured in an OMIT-style exper¬ 
iment [19j [21MTT] , where one measures the reflection of a 
weak additional probe tone incident on the cavity. The 
cavity spectrum in contrast also yields information on 
polariton occupancies. As the polariton energies E„ are 
well separated, the output spectrum will have a series of 
peaks corresponding to the emission or absorption of a 
given polariton species. The magnitude of these peaks is 
directly proportional to the occupancy of the given po¬ 
lariton. 

It is also useful to look at the total number of pho¬ 
tons due to a given polariton resonance, which we can 
obtain by integrating the output spectrum around the 
corresponding resonance. We thus introduce: 

i-Uo+Su i 

n^[w 0 M= (63) 

where ujq will be taken to be E± , and Sui will be taken to 
be larger than the spectral width of the given polariton 
resonance. 
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In what follows, we consider signatures of quan¬ 
tum heating in the spectrum for particularly interesting 
choices of the drive laser detuning A. 


B. Limit of zero temperature dissipation: effects of 
nonlinear interaction on quantum heating 

The first studied limit is for a mechanical bath at zero 
temperature ( n = 0), where the finite number of polari- 
tons inside the optomechanical cavity exclusively comes 
from quantum heating. We focus on two regimes: the 
red sideband drive, i.e. A = , where both polari- 

tons are equal mixture of photons and phonons, and the 
asymmetric polaritons regime, where the polaritons are 
not equal combinations of photon and phonon. For the 
latter regime, we chose A = — 1.8u >m as a representa¬ 
tive laser detuning. The results predicted for the cavity 
driven on the red sideband has the advantage to be ro¬ 
bust to temperatures since both polaritons have an im¬ 
portant photon part; this implies that even for finite 
quantum heating is still the prevalent source of polari¬ 
tons. In contrast, a laser detuned at A = —1.8leads 
to more striking modifications of the polaritons energy 
distribution since C° ff [A = —1.8 com] > C® ff [A = -wn] 
(cf. Fig. [5]). In both cases, we show that the nonlinear 
interaction modifies the energy distribution of the polari¬ 
tons as it tends to thermalize the two species. 


1. Results for symmetric polaritons (red sideband drive) 


In Fig. [8] we plot the cavity DOS (Eq. (|60|), the cav¬ 
ity spectrum (Eq. (59)), its energy distribution function 


(Eq. rt6Tj) ) and the corresponding effective temperature 
(Eq. ( |62| ) for A = —ojm- The linearized theory (g = 0) 
is compared to the case where g = n. For the latter 
interacting case, we present results obtained from three 
different methods: the leading order in g self-energies 
(Eqs. (42) and (52)), the self-consistent approach de¬ 


scribed in Sec. IV B and finally, a numerical simulation 
of the Lindblad master equation given in Eq. (21). By 


comparing the three different approaches, one sees that 
for g = k, higher order corrections captured by the self- 
consistent approach play an important role for the effec¬ 
tive distribution functions and the effective temperatures 
(panel (e)-(h) of Fig. 8|. In contrast, the DOS is already 
well described at the leading order in g , which is in agree¬ 
ment with Ref. HU 

The splitting of the + polariton resonance in pd[u] and 
S'dfw] near w = E + arises from the hybridization between 
the states |+) and |—, —); the resulting hybridized states 
become spectrally resolved for 5 > k (see HU for more 
details). The same hybridization phenomena gives rise 
to a resonance in and Tf l [ui) at the — polariton 

frequency, w = E- (see Fig. |8](e) and (g)). 

In Fig. [9j we plot the number of photons leaving 
the cavity near each polariton resonances, as defined in 


Eq. (631, and show that one can effectively isolate the 
effects of nonlinear interaction. To do so, one varies G 
around G res such that the nonlinear interaction get am¬ 
plified by a factor of um/k when the nonlinear process 
becomes resonant, i.e. for G = G res , compared to the off- 
resonant case, i.e. G — G res > k. Away from resonance, 
the number of photons leaving the cavity is in good 
approximation given by the linearised theory (dashed 
lines in Fig. [9]). For the parameters here (n^ = 0 and 
7 /k < 1), the linearised theory leads to (see Eqs. (16) 
and (A4|): 


< 64 » 

Fig. [9] clearly shows the thermalization between polari¬ 
tons brought about the nonlinear interaction. Without 
interactions, at A = -% the — polaritons have a lower 
effective temperature than the + polaritons, c.f. Fig. [3](c) . 
When G is near G res , the nonlinear interaction “turns on” 
and allows the two polariton species to exchange energy 
and partially thermalize (i.e. interactions heat up the — 
polaritons while cooling down the + polaritons.) 


2. Results for asymmetric polaritons 


In Fig. 

in the case where the laser detuning is A = — 1.8o;m (and 
again, G = G res ). For this more negative detuning, the 
polaritons are no longer an equal mixture of photons and 
phonons: the - polariton is more phonon-like while the 
+ is more photon-like. This particular asymmetry leads 
to larger values of C° ff than in the red sideband regime, 
mainly because of the long-lifetime of the phonon like po¬ 
lariton (cf. Fig. [H]). Due to these larger C® ff , the results 
obtained to the leading order in g are not sufficient to 
recover the numerical simulation of the Lindblad master 
equation (cf. Eq. pT] )) even for the DOS. It is than crucial 
to use the self-consistent approach to properly describe 
the effects of nonlinear interaction. Moreover, from pan¬ 
els (g ) and (h), one sees that unlike the case A = — um, 
the nonlinear interaction cools down the - polaritons and 
heats up the + polaritons. This is also consistent with a 
partial thermalization, as for A = —1.8 com, without in¬ 
teractions the effective temperature of the — polaritons 
is greater than that of the + polaritons (c.f. Fig. [3jc). 


10 we plot the same functions as in Fig. [8j but 


C. Parametric amplification of the — polaritons 


While our emphasis in this work has been on quantum 
heating effects involving zero temperature dissipation, 
our approach can also conveniently describe thermal non¬ 
linear phenomena. Perhaps the most striking example of 
this occurs for detuning A ss -wm/ 2, where leading or¬ 
der perturbation theory predicts the presence of a para¬ 
metric instability at finite temperatures (see Sec. VD). 
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FIG. 8: (Color online) (a) and (b) Cavity DOS near the — 
polariton resonance and the + polariton resonance respec¬ 
tively for a cavity driven on the red sideband and when the 
nonlinear interaction is resonant (i.e. G = G res ). We work 
in the frame rotating at the drive frequency so that u> = 0 
refers to the drive frequency while E+ = 2 E- = 63.24k in 
this frame. The light dashed curves represent the linearized 
theory, the green dotted ones represent the results to lead¬ 
ing orders in g (Eq. (421 for the DOS), the full blue curves 
are the results of the self-consistent approach as described in 
Sec. IIV Bl and the black curves are for the numerical simula¬ 
tion of the Lindblad master equation shown in Eq. ( |21| ) . (c) 
and (d) Cavity spectrum in the same conditions; the results to 


leading orders in g are given in Eqs. (521 and (591, (e) and (f) 
Cavity energy distribution function (Eq. (6li)and, (g) and 
(f) the corresponding (frequency dependent) effective temper¬ 
atures (Eq. ©)• The parameters used for all the curves are 
7 /k = 10~ 4 and o;m/k = 50, which leads to the leading or¬ 
der effective cooperativities (i.e. Eq. (41)) Ci ff = 0.18 and 
Cf = 2.46. 


In that regime, the nonlinear interaction acts as an in¬ 
coherently pumped degenerate parametric amplifier with 
the + polariton (mainly phonon) being the (incoherent) 
pump and the — polariton (mainly photon) being the 
signal mode. 

the 


In Fig. 11 we show the cavity DOS (cf. Eq. pOj 
cavity spectrum Sd[u;] (cf. Eq. (59)) and the correspond¬ 
ing energy distribution function (cf. Eq. fl6l])) for 

a laser detuning A = —0.65 wm and «th = 650. Com¬ 
bined with a damping rate of the mechanical resonator 
7 = 10~ 3 k, its resonant frequency ojm = 50k and g = k, 
one obtains, to leading order in an effective coope rativ- 


VD 


ity C5 = —0.97 (cf. Eq. (411). As discussed in Sec. 
for CI ff = — 1, the leading order perturbation theory pre¬ 
dicts a parametric instability caused by the nonlinear in¬ 
teraction. For CO® > — 1, one thus expects an important 



FIG. 9: (Color online) (a) Output photon flux in a band¬ 
width of 5k near the - polariton resonance E_ (c.f. Eq. ( |63| ) ) 
as a function of the many photon coupling G, and for a con¬ 
trol laser detuning A = —ujm- (b) Same as (a) but near the 
+ polariton resonance E+. Note that G can be varied by 
simply tuning the amplitude of the control laser. On reso¬ 
nance, i.e. G = G r es, the nonlinear effects are enhanced by 
a factor of cum/k = 50 compared to the off-resonance case, 
i.e. (G — G res ) > k. The light dashed curves represent the 
linearized theory (g = 0), the solid blue curves are for the self- 
consistent approach as described in Sec. |IV B| and the black 
dashed curves are for the numerical simulation of the Lind¬ 
blad master equation shown (211, but this time, using the full 
nonlinear part of the Hamiltonian in Eq. © that includes all 
the non-resonant nonlinear processes. For all the curves, we 
used 7 /k = 10~ 4 and njf; = 0. 


narrowing of the cavity DOS as well as an important 
heating of the cavity near the — polariton resonance. 
This predictions from the leading-order self energy are 
shown in Fig. |TT] 

Not surprisingly, higher-order corrections (as captured 
by the self-consistent self energy) are especially impor¬ 
tant in this regime and strongly contribute to prevent 
the system from going unstable. More precisely, it is the 
hybridization between the states |+) and |-) that com¬ 

petes with the parametric amplification of the - polari- 
tons; the high number of - polaritons leads to an impor¬ 
tant modification of the energy of the hybridized states 

^(1+) ± |-)) so that the nonlinear interaction cease 

to be resonant. The result of this competition is shown 
clearly in Fig. While the leading order theory pre¬ 
dicts a photon occupancy at the — polariton resonance of 
n e ^[E_] « 50, in the self-consistent approach, one only 
obtains a value ~ 1.3. This is of course still much larger 
than what would be obtained without interaction; in that 
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FIG. 10: (Color online) Same as Fig. [ 8 ] but for laser detuning 
A = —1.8 uim- In that case, E+ = 2 E- = 92.08k in the rotat¬ 
ing frame and the corresponding leading order (i.e. Eq. (411) 
effective cooperativities Cl* = 1.92 and Cf = 5.40. 


case, nf[EJ\ ~ 0. 


D. Effective two phonons absorption 



FIG. 11: (Color online) Same as Fig. [ 8 ] but for laser detuning 
A = — 0.65wm and finite temperature = 650. In these cir¬ 
cumstances, leading order perturbation theory predicts that 
the system is close to a parametric instability (see Sec. V D l. 
Note the mechanical damping rate here is larger than in previ¬ 
ous plots (7 = 10 _3 k), as this enhances the pumping effects of 
mechanical temperature. For u>m = 50 k and g = k, the lead¬ 
ing order effective cooperativies (cf. Eq. @) are Cl* = -0.97 
and Cf = 0.75. In the frame that rotates at the drive fre¬ 
quency, E+ = 2 E- = 53.34k. One sees from these plots the 
crucial importance of higher-order-in -7 terms. 


Another striking example of temperature-enhanced 
nonlinear effects occurs for a laser detuning A « —2 u>m- 
In this regime, the — polaritons are mostly phonons, 
which leads to a thermal “stimulated emission” enhance¬ 
ment of S^[w], leading to a strong modification of the 
cavity DOS. This physics follows from Eqs. (39), and was 
discussed extensively in Ref. [T7J However, as shown in 
Fig .0 one also obtains significant nonlinearity-induced 
heating of the cavity in this regime, as we now describe. 

Recall that if A = —2 u>m, then G res = 0 and the + 
(—) polaritons are exactly photons (phonons). A neces¬ 
sary consequence is that the amplitude for the resonant 
nonlinear interaction vanishes, g = 0. One thus ideally 
wants a detuning close to, but not exactly equal to — 2 u>m 
such that 0 < G res <C ojm- In the high temperatures limit 
(h^ 1) and for weak nonlinear interaction Cf <C 1 

(i.e. g < k), the expression to lowest order in G/ojm for 
the + resonance in the cavity spectrum simplifies to: 


Sd[u] ~ - (—V (-) 

7 V UJ M / \k/ 


16 «) 


8 g^k\ 3 (w - E + ) 2 + k 2 _ ' 


(1 4. 8GiO 

V + 7j 


(65) 


As discussed near Eq. (25), the competition between 


the parametric heating associated with the linearized op¬ 
tomechanical interaction and the standard optomechan¬ 
ical optical damping leads to an optimal value of G/lum 
where the heating is maximal. In the optimal case, 
(G/w«) 2 = and ' S 'd[- B +] ~ Details of 

the calculations that leads to Eq. ( [65]) are presented in 
Appendix [Cj 

In Fig. [l2j the cavity DOS and the cavity spectrum 
near the + polariton resonance is plotted for a laser de¬ 
tuning near — 2 ojm and for finite mechanical bath tem¬ 
peratures n^. In panel (a), we show the sharp dip in 
the density of state due to nonlinear interaction, also de¬ 
scribed in Ref. H3 Note however that effects of higher or¬ 
der in 3 , captured in the self-consistent approach, consid¬ 
erably modify the predictions made in Ref. II71 where only 
effects to leading order in g were considered. This sharp 
feature is completely analogous to the optomechanical 
induced transparency (OMIT) observed in the optome¬ 
chanical cavity weekly driven on the red sideband {19[ F29V - 
EB, except that here, it is the nonlinear interaction that 
is involved. In panel (b), we show that nonlinear interac¬ 
tion greatly modifies the cavity spectrum; no output light 
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c j/k 


FIG. 12: (Color online) (a) Signatures of the effective two- 
phonon absorption in the cavity DOS near the + resonance; 
the inset shows a larger range of frequency than the main plot. 
The light-blue dashed curve represents the linearized theory 
(g = 0 ), the gree n do tt ed o ne rep rese nts the results to leading 
order in g (Eqs. (421, (521 and (59») and, the full dark-blue 
and black curv es are the results of the self-consistent approach 
(S.C.) of Sec. (IV B) for = 100 and = 0 respectively, 
(b) Signatures of the effective two-phonon absorption in the 
cavity spectrum in the same circumstances than (a). It shows 
the striking temperature-enhanced effect of nonlinear interac¬ 
tion in the + polariton population. The parameters used are 
■y/n = 10 -4 and uim/k = 50. 


is predicted by the linearized theory while a sharp signal 
is produced when nonlinear interaction is considered. 

The effective two phonons absorption becomes a very 
important process in the high temperature limit. In that 
case, the nonlinear interaction becomes easier to observe, 
but the phenomenon tends to become purely classical. In 
order to support this statement, we present a classical 
treatment of this phenomenon in Appendix [D] and show 
that it succeeds to recover the right dependence in tem¬ 
peratures and single-photon coupling constant g of the 
two phonons absorption signature in the cavity spectrum 
at the lowest order in g. 

VII. CONCLUSIONS 

In this work, we have described the effects of nonlin¬ 
ear interaction on the non-equilibrium state of the op¬ 
tomechanical cavity. We have shown the tendency of the 
nonlinear interaction to thermalize the polaritons, that it 
can lead to a new parametric instability for a red-detuned 


laser drive as well as a temperature-enhanced effective 
two phonons absorption. In addition to these results, 
we have presented in details many technical aspects with 
the aim to provide the proper tools to investigate nonlin¬ 
ear effects in more complicated optomechanical systems. 
This work also opens the path to a more detailed charac¬ 
terization of this new parametric instability and further 
investigations of its consequences. 


Acknowledgements 

We thank Nicolas Didier and Felix Beaudoin for help¬ 
ful discussions. This work was supported by NSERC and 
FQRNT. 


Appendix A: Change of basis for a detuning on the 
red sideband (A = — ujm ) 

In this Appendix, we show particular examples of the 
transformation used to go from the photon/phonon basis 
to the polariton basis for A = — ujm■ 

We start with the change-of-basis coefficients a b / d ,± 
and a b / dj ± introduced in Eqs. (|5j), 


1 UJM + E± 

a b ,± = ± ~ - r— , 

V™M y/E± 

(Al) 

, 1 wm - E± 

/q / y—, 5 

V»WM y/E± 

(A2) 

1 UJM + E± 

(A3) 

/Q- pp I— 5 

V°UJM y/E± 

1 wm - E± 

/q / i 

V°UJ M y/E± 

(A4) 

with 


E± = ujm\J 1 ± 2 G/ujmi 

(AS) 

Gres/WM = 3/10. 

(A6) 

From these coefficients, we can express the different ef- 

fective nonlinear coupling constants g aa ^ a n 

and A± in- 


troduced in Eq. ([6]). Here, we only show few examples. 

<#+_+$+-++A+ 


=g [a dt+ a dt+ (a b - + a b _) 

+ (dd,-CTd,+ + &d,+0!d,-)(o? b: + + db,+ )] , 


/ \ 3 / 4 
g (um\ 

' E_ 

- E + + ^( E -+‘ 2E + )] 

fcal 

i 

E+ 


(A7) 


g = 9 -—f — 9 [ a d-(^d,-{o/-b,+ + cq>,+ ) 

+(ad,-&d,+ + otd,-Ud,+ ){ot b ,- + <W_)], 


9 

^m\ 3/4 

\ (1 + 2+ (1 - - g -)E +1 

4^2 

K E +) 

£L 

(AS) 
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A- —g [( 2 otd,-&.d.- + cy^ _ + + &b,— ) 

+ {6td.,-Otd.,+ + &d,+Oid,-){o!b,+ + £*&, + )] ) 


9 

fcj M \ 3/4 

r2(— - 1 )E+ + -^-(E+- £L) 1 

V CO M ' ^ \ T- / 

4^2 

KE-) 

1 

+ 


(A9) 


Appendix B: Parametric amplifier description 


In this appendix, we derive the condition to get an in¬ 
stability from the parametric amplifier Hamiltonian and 
get the corresponding mean number of signal excitations. 
These calculations are based on the formalism introduced 
in [50] . 

Starting with the nonlinear interaction in the effective 
Hamiltonian of Eq. ([ 8 ]) and considering the hypothetical 
case where the + mode is coherently pumped, such that 
(c_|_) = n° + , we get, in the mean field approximation: 


along this work, we find that the mean number of — 
polaritons inside the optomechanical cavity is 


(c_c_) 


\Q\ 2 

i-IQI 2 


1 6g 2 n+ 
k 2 _ — 16 g 2 n° 


(B8) 


From this result, one immediately sees that for \Q\ 2 = 1, 
the mean number of polaritons diverges. This condition 
sets the threshold above which the parametric amplifier 
goes unstable. More precisely, the system becomes un¬ 
stable when 


n° + > 


16 g 


2 ' 


(B9) 


These results have exactly the same form than Eqs. (56) 


and (57), but here, we have explicitly used the fact that 
the + polaritons are coherently pumped such that (c+) = 


n+. 


H oS 


C-C- + H.c.). 


(Bl) 


Using standard input-output theory, one derives the fol¬ 
lowing equation of motion 


d t c- 


i[H eS: C-} - C— - 

(B2) 


(B3) 


Here, £_ represents an incoming field, which includes 
noise coming from the bath coupled to the — polaritons 
(as in Eq. (33)). Since we will focus only on quantities 
evaluated at ui = E_ and that H e s of Eq. (Bl) is written 
in the interaction picture where E _ = 0 , we can seek for 
the particular solution given by cLc_ = 0. In that case, 
we have: 



=* c_ = -j= (1 - IQI 2 ) 1 (q£ - £_), 


with 


Q 


4ig 



K- 


(B4) 

(B5) 


(B 6 ) 


Imposing the right commutation relation [c_,c_] = 1, 
one finds 




K - (i - IQI 2 ) 

4 


(B7) 


which has the right units since (£L£_) represents a rate at 
which excitations are coming in. In the case of Gaussian 
noise with zero mean (i.e. (£_£_) = 0 ), as studied all 


Appendix C: Expansion of the photon spectrum 
function for A S 3 —2 ujm 


In this Appendix, we show details of the calcula¬ 


tion that leads to Eq. (65) starting from Eqs. (52) and 


Eq. (59). Here, the limit of interest is for A ss —2cj m , 


G/ujm ~ j/k -C 1, 1 and C+ ff <C 1. In these cir¬ 

cumstances, photons are mainly + polaritons, such that 
Eq. (65) for frequencies near E + reduces to 


Sd[u]\u>^E + ~ 2nn f f[u]p + [ui}. 


(Cl) 


Also, for A ss —2u>m and G/ojm ~ j/k <C 1, the mean 
number of photons inside the cavity without nonlinear 
interaction, , is negligible compared to the nonlinear 


contribution. Thus, from Eq. (52), we see that n+fw] 


reduces to a sharply peaked Lorenztian with witdh 7 + 


7(1 + r opt ) ~ 7 with the optical damping rate 
, 8G 2 


K- 

r opt ~ |^-k. Consequently, p+[oS\ (Eq. (42|), which 
has a characteristic width of S3 k. can be evaluates at 
E + . Doing so, one gets 


S , d[w]| 


CU^Ea. 


2tt nf [uj\p+[E + ], 


P+[E+] = - 


n K+ (l+Cf) 7 TK 


(C2) 

(C3) 


Now, we can write as (we have dropped and 

explicitely used E + = 2E_) 


hf[u] S3 n^Cf 


(id — E .|_) 2 + K_ ' 


(C4) 


with 


= 


eff _ 4 5 2 (n °_) 2 




(C5) 
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Thus 


cm ^16 5 2 ( n °_) 2 

- k 2 k _ (w _ E+ y + f~2- > 


(C 6 ) 


Again, for this particular limit, we can approximate n_ ? 
and g 2 ~ (^yjj) 2 # 2 ) such that we get Eq. (651, i.e. 


K*Sd[w]| 




G_ 
7 V w m 


(f) ; 


16(S“) 


(1 + ^ 2 Ei) (u - E + ) 2 + k 2 _ 


(C7) 


Appendix D: Classical treatment of the effective two 
phonons absorption 


As shown in the Sec. VI D[ the nonlinear process corre¬ 
sponding to the absorption of two phonons by a photon 
from the classical drive is greatly enhanced by tempera¬ 
ture. In particular, even if the nonlinear interaction di¬ 
rectly comes from single photon dynamics, one can show 
that in the high temperatures limit, this phenomenon 
becomes purely classical. In what follows, we present an 
accurate classical description of this limit. 

The classical Langevin equations for the cavity field 
coupled to a mechanical resonator via radiation pressure 
force in the limit of weak displacements are 


ait) = - iu c ^ +ia in e lWLt , 

(Dl) 

x(t) + ^x + ujfX = . (D2) 

2 ttiL/q 

Here, m is the mass of the mechanical oscillator, L 0 is 
the length of the cavity when x{t) = 0 and the damping 
rate k has been introduced following the standard input- 
output formalism |5f)j . The cavity field is normalized 
such that the steady state mean energy inside the cavity 
is H cav = \a\ 2 = K|ai„| 2 /((f ) 2 + A 2 ), where a is defined 
as the steady state mean value of the cavity held, a[t) = 
(a + 5a{t))e~' lulLt . 

We solve these nonlinear coupled equations of motion 
via a perturbation approach. Starting with a sinusoidal 
displacement x(t) = rrosinwe get the consequent 
evolution of the optical mode and then, the perturbed 
displacement of the mechanical resonator. This is a good 
approximation only in the weak coupling regime (G -C 
k) where the eigenstates have a well defined number of 
photons and phonons. This method is thus restricted 
to the extreme detuning case (A ss —2 u>m) studied in 
Sec. Ena 


Following this method, one finds that the steady state 
cavity held is given by 


i{t) = 


+oo 

E 


n,m =—oo 


§ + i(nujM - A) 




(D3) 

Here, J n {x) are the Bessel functions of the hrst kind and 
where, in the small displacement limit, e = plays 

the role of the perturbation parameter. 

Note that by computing the radiation pressure force 
acting on the mechanical resonator in Eq. (D2) and by 
expanding to the lowest order in e, one recovers the stan¬ 
dard optical damping and optical spring constant coming 
from the linear interaction (cf. Eq. ( |2b| )) in the weak cou¬ 
pling regime |40] . 

We focus here on the peak in the cavity spectrum 
that comes from the nonlinear interaction. According 
to Eq. (58), the proper definition of the classical coun¬ 
terpart of the photon spectrum to use is 


SdM = 



dte lalt (Sa* (t 0 )Sa(t 0 + t)) to 


(D4) 


Here, the fluctuations of the cavity held S a(t) are given by 
keeping only the n ^ —to terms in Eq. ( |D3 ). The mean 
value (...)t 0 means taking the average over all initial time 
to and the coefficient 1 /up is there to ensure that units 
match the definition of the photon spectrum used in the 
quantum treatment (cf. Eq. (58)). 

Using Eq. (D3) and focusing only on the frequency 
range near 2 ojm with A ss one hnds that in the 

good cavity limit (k <C u>m) > the classical cavity spec¬ 
trum reads 


Sd[w] ~ 16 ( —) (-) (n“) 2 7r<5(w - 2 w m )- (D5) 

\WM/ ' KJ 


Where we have used the dehnition of the single photon 


coupling constant, g = 


; = IZPF LJC 


, with Xzpf being the zero 


point motion of the mechanical resonator and the semi- 
classical relation 2 = ujfnfh valid for 3> 1 . 

If one takes T opt = 0 and 7 —> 0 in Eq. (|65|), S^oj ss 
E + \ and S c i[u; ~ 2 wm] agree. The classical approach fails 
in getting the width of the peak since the mechanical 
damping rate is not taken into account while calculating 
the effects of the mechanical motion on the cavity field. 
Moreover, the optical damping doesn’t come out in this 
result in the level of approximation we have used. 
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